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oc : 

Based on the cranking model and the random phase approximation, we point out that 
the wobbling excitation on top of the s band in 182 Os is stable against angular momentum 
£SJ ' tilting. This is consistent with the general trend that the wobbling excitations in 7 < 

rotating nuclei are more stable than those in 7 > ones found in our previous studies. 
In higher N isotopes known to be 7 soft, however, a different type of tilting instability is 
expected. Its possible correspondence to the experimental data is also discussed. 

t> : 

Symmetry breaking in the nuclear mean field is an analog of second order phase 
transitions in infinite systems and is one of the key concepts in the theory of nuclear 
collective motion. According to the general concept of symmetry breaking, when 
one approaches the transition point from the symmetric side, softening of collective 
vibrational mode takes place as a precursor of the phase transition. Examples are: 
1) In the spherical to axial shape transition, the 2 + quadrupole vibration softens, 2) 
in the axial to triaxial shape transition, the 7 vibration softens, and 3) in the normal 
■ fluid to superfluid transition, the pair transfer cross section increases. 

1 ■ Collective rotation of axially symmetric nuclei takes place only about a principal 

I/-) I axis (usually named the x axis) perpendicular to the symmetry axis (the z axis). 

In triaxially deformed nuclei, however, rotations about all three principal axes are 
possible. Therefore, if triaxiality sets in gradually, the angular momentum vector 
starts to wobble when seen from the principal axis frame. Eventually the angular 
momentum vector tilts permanently from the x axis. This regime is called the tilted 
axis rotation (TAR) in contrast to the usual principal axis rotation (PAR). Thus, 
J> I the softening of the wobbling motion is the precursor of symmetry breaking from 

the PAR to the TAR. We call this instability of the PAR mean field, caused by 
the softening of the wobbling motion, the tilting instability. After this instability, 
a TAR mean field, in which the signature quantum number that is associated with 
a 7r rotation about the x axis is broken, replaces the PAR mean field. As shown in 
Eq.{H), the excitation energy of the wobbling motion is determined by moments of 
inertia, which are dynamical response of the system to rotation from a microscopic 
viewpoint. Therefore, not only moments of inertia depend on 7 deformation but also 
7 deformation itself depends on the rotation frequency since rotational alignments 
of quasiparticles exert shape driving effects on the whole system according to their 
positions in the shell. 

The small amplitude wobbling motion at high spins was first discussed by Bohr 
and Mottelson 1 ) in terms of a macroscopic rotor model with constant moments of 
inertia. Then it was studied microscopically by Janssen and Mikhailov 2 ** and Mar- 
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shalek 3 ) in terms of the random phase approximation (RPA) that gives dynamical 
moments of inertia. Since the small amplitude wobbling mode has the same quan- 
tum number, parity ir = + and signature a = 1, as the odd-spin member of the 
7 vibrational band, Mikhailov and Janssen 4 ) anticipated that it would appear as 
a high-spin continuation of the odd-spin 7 band. But it has not been clear that 
in which nuclei, at what spins, and with what shapes it would appear. Using the 
RPA, Shimizu and Matsuyanagi 5 ) studied Er isotopes with small | — y | , Matsuzaki 6 ) 
and Shimizu and Matsuzaki 7 ) studied 182 Os with a rather large negative 7 but their 
correspondence to experimental data was not very clear. In 2001, 0degard et a/. 8 ) 
found an excited triaxial superdeformed (TSD) band in 163 Lu and identified it firmly 
as a wobbling band by comparing the observed and theoretical interband E2 transi- 
tion rates. These data were investigated in terms of a particle-rotor model (PRM) by 
Hamamoto 9 ) and in terms of the RPA by Matsuzaki et a/. 10 -* In the latter, the calcu- 
lated dynamical moments of inertia are rotation frequency dependent even when the 
shape of the mean field is fixed. This dependence is essential for understanding the 
observed behavior of the excitation energy. In 2002, two-phonon wobbling excita- 
tions were also observed by Jensen et aZ. 11 ) and their excitation energies show some 
anharmonicity. In Ref. 10), a numerical example of the softening of the wobbling 
motion in the positive-7 nucleus, 147 Gd, was presented. Matsuzaki and Ohtsubo 12 ) 
elucidated that study by examining shape change of the potential surface as a func- 
tion of the tilting angles. In that paper it was also discussed that the observed 
anharmonicity may be a signature of the onset of softening. Oi 13 ) proposed a new 
model to account for this softening. Almehed et a/. 14 ) also discussed this. Recently 
Tanabe and Sugawara-Tanabe proposed an approximation method to solve the PRM 
and applied it to the TSD bands. 15 ) Kvasil and Nazmitdinov gave a prediction for 
the wobbling excitations in normal deformed nuclei 16 ) by utilizing the sum rule type 
criterion found in Ref. 17). 

The excitation energy of the wobbling motion is given, as a function of moments 
of inertia, by 1 ) 

hw woh = huj rot > y — , (l) 

y UyU z 

where uj TO t is the frequency of the main rotation about the x axis and Js are moments 
of inertia about three principal axes. This indicates that J x > Jy, Jz or J x < J y , J z 
must be fulfilled for co wo b to be real. The irrotational model moment of inertia is 
given by 

j; r KS in 2 ( 7 + ^), (2) 

with k = 1 — 3 denoting the x — z principal axes, and its 7 dependence is believed to 
be realistic. When this is taken, —60° < 7 < for the former or —120° < 7 < —90° 
or 30° < 7 < 60° for the latter is required. Since the 7 deformation of the observed 
TSD band is 7 ~ +20°, another mechanism is necessary for the wobbling excitation 
to exist. It was found in Ref. 10) and elucidated in Ref. 17) that the alignment of the 
last odd quasiproton brings an additional contribution to J x and consequently makes 
Jx > Jy m place of J x < J y in the irrotational-like behavior. But the smallness of 
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J x — Jy implies fragility of the excitation. 

The negative-7 collective rotation, —60° < 7 < 0, is expected to occur prevail- 
ingly. However, it looks difficult to excite a wobbling mode on the ground band 
of even-even nuclei because of J x ~ J y in those cases (see Ref. 6)). Therefore, 
the following three conditions are desirable for the wobbling excitation to exist: 1) 
—60° < 7 < 0, 2) I7I is not small, and 3) existence of aligned quasiparticle(s) that 
makes J x larger. Prom these conditions, we chose the s band of 182 0s as a repre- 
sentative in Refs. 6), 7). We concluded that a wobbling excitation exists on top of 
the s band of 182 Os. Recently, Hashimoto and Horibata 18 ) presented the opposite 
conclusion. Here we briefly comment on their work before proceeding to the main 
discussion of this paper. They recently reported a renewed three-dimensional crank- 
ing calculation for 182 Os paying attention to the stability of the s band based on 
their previous calculation. 19 ) They concluded that the wobbling excitation on the s 
band does not exist; this contradicts our previous calculation. 6 )' 7 ) A close looking 
into their works leads one to find that the character of the s band is different between 
theirs and ours. Although not stated in Ref. 18), it was reported in Ref. 19) that their 
s band consists of two aligned quasiprotons. Their low- J? hg/2 character would lead 
to a positive-7 shape. Note that their convention for the sign of 7 is opposite to the 
Lund convention adopted here. As stated above, wobbling excitations in positive-7 
nuclei are fragile. Although our calculation adopted fixed mean field parameters, we 
conformed to the experimental information that suggests the s band consists of two 
aligned £13/2 quasineutrons ?°> Since the Fermi surface is located at a high position 
in the i 13 / 2 shell, the alignment leads to a negative-7 shape. As discussed above, the 
wobbling excitations on negative-7 quasiparticle aligned configurations are rather 
stable. This is the reason why the conclusions of Hashimoto and Horibata and ours 
are different. The collective excitation on the g band is expected or exists in both 
calculations, but in our calculation it is 7 vibration-like rather than wobbling-like 
(see Ref. 4)). 

Now we proceed to present a numerical example of another type of tilting insta- 
bility in negative-7 rotating nuclei, different from that in positive-7 cases discussed in 
our previous works, 10 )' 12 ) although wobbling excitations are stable in many negative- 
7 cases when it exists. The meaning of "different" is elucidated later. First, we review 
our model briefly. We begin with a one-body Hamiltonian in the rotating frame, 



hmi = ^ + ^M(u; 2 x x 2 + u 2 y 2 + u 2 z z 2 ) + v ls \ ■ s + ^(l 2 - (1 2 )jv oso ). (5) 



In Eq.Q, t = 1 and 2 stand for neutron and proton, respectively, and chemical 
potentials A r are determined so as to give correct average particle numbers (N T ). The 
oscillator frequencies in Eq.([5]) are related to the quadrupole deformation parameters 
€2 and 7 in the usual way. They are treated as parameters as well as pairing gaps 
A T . The orbital angular momentum 1 in Eq.(|5j) is defined in the singly stretched 

coordinates x' k = \ h^Xk and the corresponding momenta, with k = 1 — 3 denoting 
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x - z. Since h' conserves parity tt and signature a, nuclear states can be labeled by 
them. We perform the RPA to the residual pairing plus doubly stretched quadrupole- 
quadrupole (Q" ■ Q") interaction between quasiparticles. Since we are interested in 
the wobbling motion that has a definite signature quantum number, a = 1, only two 
components out of five of the Q" ■ Q" interaction are relevant. They are given by 



(-) 

int 



K=l,2 



K ' 



where the doubly stretched quadrupole operators are defined by 

Q'k = Qk{xu -»■ x'l = —x k ), 

UJQ 



and those with good signature are 



V2(l + 6ko) 



(Qk ± Q. K ) . 



(6) 



(7) 



(8) 



The residual pairing interaction does not contribute because P T is an operator with 
a = 0. The equation of motion, 



h' + H { ~ ) X ] 



RPA 



hw n xl, 



for the eigenmode 



(a=±l/2) 



leads to a pair of coupled equations for the transition amplitudes 
Then, by assuming 7 / 0, this can be cast 3 ) into the form 



o ( ~ ] x f 



(9) 



(10) 



(11) 



2 2 



rot" 



■Jy \^n)Jz \^>r. 



= 0. 



(12) 



This expression proves that the spurious mode {u n = uj ro t', not a real intrinsic ex- 
citation but a rotation as a whole) given by the first factor and all normal modes 
given by the second are decoupled from each other. Here J x = {J x )/u TO t as usual 
and the detailed expressions of Jy^P (to n ) are given in Refs. 3), 6), 7). Among normal 
modes, one obtains 



^wob — <^rot 



Jx Jy ^(^wob)^ {^Jx Jz ^(^Nvob) 

Jy^ (w wo b) l 7i Cff ' ) (Wwob) 



(13) 
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by putting uj n = w wo b- Note that this gives a real excitation only when the argument 
of the square root is positive and it is non-trivial whether a collective solution appears 
or not. Evidently this coincides with the form (JTJ) derived by Bohr and Mottelson in a 
rotor model 1 ) and known in classical mechanics. 21 ) Further, this makes it possible to 
describe the mechanism of the tilting instability in terms of the dynamical moments 
of inertia. The wobbling angles that measure the amplitude of the vibrational motion 
of the angular momentum vector around the x axis are defined by 



i yVy^c^wob)! 2 + i-/j PA ^wob)i^ 

(i PA) ) 

Jz \^wob) 



?wob = tan -i ^ , (14) 

(Jx ) 



¥>wob = tan 1 



Jy \ui wo b) 



(15) 



with (PA) denoting the principal axis frame. The PA components of the angular 
momentum vector are defined by 

(J X PA) ) = (J X ), (16) 
iJ y FA) = tJ v ~ J (1) Q2 \ ( 17 ) 

4 PA) =Jz- ri ,X x) - {+) Q[-\ (18) 

in terms of the RPA matrix elements of their uniformly rotating frame components 
usually calculated in the cranking model, 3 )' 6 )' 7 ) because the PA frame is determined 
by diagonalizing the quadrupole tensor Q^ PA ). 3 )' 22 ) 

We choose 186 0s, bearing possible correspondence to the experimental data in 
mind. The s band consists of (i^i3/2) 2 - In this calculation, we concentrate on the 
direct rotational effect by ignoring the effect of the possible rotational shape change. 
The adopted mean field parameters are ti = 0.205, 7 = —32°, and A n = A p = 0.4 
MeV. Calculations are performed in the model space of five major shells; N osc = 3 — 
7 for neutrons and 2-6 for protons. The strengths of the 1 • s and l 2 potentials are 
taken from Ref. 23). Figure [Q reports the excitation energy fko wo ^ in the rotating 
frame. Decrease of this quantity signals the instability of the principal axis rotating 
s band that supports the small amplitude wobbling excitation. Figure [2] graphs the 
wobbling angles # wo b an d Vwob • While the angular momentum vector wobbles around 
the x axis with # wo b ^ 15° up to just below the instability point, ip wo b increases 
gradually. This means that the z component increases gradually. Eventually at the 
instability point the angles look to reach # wo b > 45° and (p wo b = 90°, that is, the 
angular momentum vector tilts to the x — z plane. Although the present calculation 
can not go beyond the instability point, a numerical example of the correspondence 
between the instability of the PAR and the TAR that follows it was presented in 
Ref. 12). More direct information about the shape that the system would favor can 
be obtained from the moments of inertia shown in Fig. [3l This figure shows that 
Jx = Jz is realized at the instability point; this is a different type of tilting instability 
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from that observed in 7 > nuclei that is caused by J x = J y . Here we elucidate the 
meaning of "different type". The instability brought about by J x = J y discussed in 
Refs. 10), 12) and that by J x = J z discussed here are similar in the sense that the 
energy costs of rotations about two different axes coincide. But here we base our 
discussion on the physical picture that 7 > and 7 < are different rotation schemes 
and, according to the reason discussed above, at 7 > 0, w wo b can not be real without 
aligned quasiparticle that makes J x larger, in contrast, at 7 < 0, <^ wo b can be real 
without it. Note that nothing peculiar happens at J y = J z because the instability 
is given by zeros of Eq. (|13p . Although selfconsistent shape change is beyond the 
scope of the present simple-minded calculation, J x = J z may indicate that either a 
TAR {J y 7^ 0) or another PAR, that is, an oblate collective rotation {J y = for the 
irrotational rotor), would be favored. A possibility of oblate collective rotation was 
first discussed by Hilton and Mang 24 ) for 180 Hf, and very recently by Walker and 
Xu 25 ) and Sun et al. 26 ) for 190 W. In the present case, J y is decreasing but not 0. 
Therefore, it is natural to regard the rotation scheme just after the instability as a 
TAR. 

0.5 J , , , , 1 

0.4 
1 0.3 




0.35 



ha> l0t (MeV) 

Fig. 1. Rotation frequency dependence of the wobbling excitation energy on the s band of 186 Os. 

Although the quantitative criterion for the occurrence of the instability is be- 
yond the scope of the present calculation, we confirmed that the instability occurs 
at lower rotation frequency for smaller £2 or larger N. These results point to con- 
sistency with the N dependence of the 7 softness in this mass region seen in the 
quadrupole deformation, 27 ) the excitation energy of the 7 vibration, ^ and the high- 
K isomerism. 28 ^ 

Finally we mention possible correspondence to the observed data. In Ref. 29), 
Balabanski et al. reported an anomalous termination of the yrast band of 186 0s 
at 18 + . According to their calculation, the {vivz/2) 2 alignment drives the shape 
to 7 ~ —30° before this termination. Actually the mean field parameters of the 
present calculation were chosen conforming to this. As for the termination itself, 
they discussed using a total Routhian surface calculation that it is related to a 
further shape change in the 7 direction. Later Wheldon et al. 28 ' discussed that it 
does not terminate. Aside from the different conclusions about the fate of the higher 
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Fig. 2. Rotation frequency dependence of the wobbling angles. 




1 1 1 1 1 1 

0.1 0.15 0.2 0.25 0.3 0.35 

hcorot (MeV) 

Fig. 3. Rotation frequency dependence of the moments of inertia. 

spin states, the yrast band changes its character at 14 + in both studies. Wheldon et 
al. 28 ^ concluded that this is caused by the crossing with the 10 + band that is tilted 
("t band"). Since the main component of the high spin part of the ground state (gs) 
band is thought to be a PAR triaxial (^13/2) 2 s band, 30 ) the observed crossing is 
ascribed to the instability of the PAR mean field qualitatively. On the other hand, 
since the 14 + and the 12 + members of the gs band correspond to fkJ ro t = 0.389 
MeV and 0.356 MeV, respectively, the observed crossing takes place between them. 
Therefore, quantitative correspondence with the present calculation in which it takes 
place at around huj TO t = 0.310 MeV is insufficient. 

To summarize, in this paper, first we have pointed out that the wobbling exci- 
tations on 7 < quasiparticle aligned bands are expected to be more stable than 
those on 7 > ones as found in our previous studies. In relation to this, we have 
clarified the reason for the different conclusions about the existence of the wobbling 
excitation on top of the s band of 182 Os between the recent work of Hashimoto and 
Horibata 18 ) and ours. Second, we have discussed, in spite of this, that the wobbling 



8 



M. Matsuzaki 



excitation in 7 < nuclei can become unstable by presenting a numerical example, 
although the quantitative criterion for the occurrence of this type of tilting instability 
is deferred to more elaborate calculations. Possible correspondence of this example 
to the experimental data is also discussed. 
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